Heat transport in disordered quantum harmonic oscillator chains 
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We study heat conduction in quantum disordered harmonic chains connected to general heat reser- 
voirs which are modeled as infinite collection of oscillators. Formal exact expressions for the thermal 
current are obtained and it is shown that, in some special cases, they reduce to Landauer-like forms. 
The asymptotic system size dependence of the current is analysed and is found to be similar to the 
classical case. It is a power law dependence and the power depends on the spectral properties of the 
reservoirs. 
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In an earlier paper [l]we studied the problem of heat 
conduction in a classical disordered harmonic chain that 
is connected to heat baths specified by general spectral 
functions. In this paper we use the formalism, developed 
by Ford, Kac, and Mazur [2] (FKM) to model quantum 
mechanical heat reservoirs, and extend our earlier calcu- 
lations to the quantum mechanical regime. Our results 
should be useful in interpreting recent experiments [3] on 
heat transport in insulating nanowires and nanotubes. 
They are also of interest in the context of the question 
of validity of Fourier's law in one-dimensional systems, a 
problem that has received much attention recently [4] . 

As pointed out in an earlier paper [9], on transport in 
Fermionic wires, the FKM method seems to be a pow- 
erful and relatively straightforward method for studying 
nonequilirium transport. Here we use this formalism to 
obtain exact results for the heat current in a disordered 
harmonic chain. Landauer's result, suitably generalized 
to the phononic case, is obtained as a special case. For 
general reservoirs the conductance of a wire depends on 
details of the reservoirs and contacts. Our analysis here 
closely follows the one followed in [1,5,6] for the case of 
classical oscillator chains. The results obtained in the 
classical oscillator case and the present problem have a 
lot of formal similarity which of course is because both 
the problems are linear. There has been some earlier 
work on quantum wires [7,8] which follow a similar ap- 
proach but we give a more clear and complete picture 
and make some interesting predictions for experiments. 

We consider a mass disordered harmonic chain con- 
taining N particles with the following Hamiltonian: 
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H = 



i=i 



N-l 



{xi -x l+1 ) 2 (x 2 +x 2 N ) 
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where {xi} and {pi} are the displacement and momen- 
tum operators of the particles and {mi} are the random 
masses. Sites 1 and N are connected to two heat reser- 
voirs (L and R) which we now specify. We model each 
reservoir by a collection of M oscillators. Thus the left 
reservoir has the following Hamiltonian: 
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where Ki m is a general symmetric matrix for the spring 
couplings, {Xi,Pi} are the bath operators and {Xi,Pi} 
are the corresponding normal mode operators. They are 
related by the transformation Xi = ^ s Ui s X s where 
Uis, chosen to be real, satisfies the eigenvalue equa- 
tion J2i K mUi s = Lo 2 s U ns for s = 1, 2...M. The an- 
nihilation and creation operators a s , a\ are given by 
a s = (P s — ilo s X s ) I '(^ojs) 1 / 2 , etc. and n s — a\a s is the 
number operator. 

The two reservoirs are initially in thermal equilibrium 
at temperatures Tl and Tr. At time t = r the system, 
which is in an arbitrary initial state is connected to the 
reservoirs. We consider the case where site 1 on the sys- 
tem is connected to X p on the left reservoir while N is 
connected to X p > on the right reservoir. Thereafter the 
whole system evolves through the combined Hamiltonian: 

Ht = H + Hl + Hr — kxiXp — k'xjsiXpi . 

The Heisenberg equations of motion of the system vari- 
ables are the following (for t > r): 



in 



\X\ = -[2xi 



X2\ 



kX ri 



mm = — (— ij-! + 2xi - xi + i) 1 < I < N 
m N XN = — [— scjv-i + 2x N ] + k'Xpi. (3) 



We note that they involve the bath variables X p y . How- 
ever these can be eliminated and replaced by effective 
noise and dissipative terms, by using the equations of 
motion of the bath variables. Consider the equation of 
motion of the left bath variables. They have the form: 



X n 

X n 



-K nl Xi n^p 
-KpiXi + kxi 



(4) 



This is a linear inhomogeneous set of equations with the 
solution 



1 



X n = - t)Xi(t) + G nl (t - t)Xi(t) 

i 

/oo 
dt'G np (t-t')kxi(t') where (5) 

F nl (t) = 6(t)^U ns Ul a 
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Thus wc find that X p (say) appearing in Eq. (3) has 
the form X p (t) = h(t) + k dt'G pp (t - t')ari(f')- The 
first part, given by h(t) = X)jF p/ (t - t)Xi(t) + G pl (t - 
t)Xj(t), is like a noise term while the second part is like 
dissipation. The noise statistics is easily obtained using 
the fact that at time t = r the bath is in thermal equi- 
librium and the normal modes satisfy (a|(r)a s /(r)) = 
f(u> s , (3l)S S s' ■ Here / = l/(e' 3 " — 1) is the equilibrium 
phonon distribution and (O) = Tr[pO] where p is the 
reservoir density matrix and TV is over the reservoir 
degrees of freedom. We define the Fourier transforms: 
= ^J^dtxtW*, G+» = JZedtG^ty-t 
and h(u>) = ± dtfi(t)e iwt . Taking limits M — ► oo 
and t — ► — oo wc get: 



X p (u) = h{uj) + kG+ p {uj)xi{uj) 
{h(u;)h(uj r )) = I(u)5(u + uj') where 



(6) 
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b(w) = V" — — - lu s ) - S(cj + u) s 
— ' loj- 



Similarly for the right reservoir we get X p > = h'(u>) + 
k'Gp, p ,(uj)xN(uj), the noise statistics of h'(u)) being now 
determined by /?'. The left and right reservoirs are in- 
dependent so that {h(oj)h' (u 1 )) = 0. Wc can now obtain 
the particular solution of Eq. (3) by taking Fourier trans- 
forms and plugging in the forms of h(u>) and h'(u>). We 
then get: 



iuit 



(7) 
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Z = 4> lm - Ai m with 

4>lm = -(8l,m+l + 8l,m-l) + ( 2 ~ ™i^ 2 )<Vr, 

A lm = 6i,m[k 2 G+ p (Lj)6i tl + k l2 G+ pl (oj)5 UN ] 
hi(co) = kh(w)S^i + k'h' (uj)8^n- 



We can now proceed to calculate steady state values of 
observables of interest such as the heat current and the 
temperature profile. We first need to find the appropri- 
ate operators corresponding to these. To find the cur- 
rent operator j we first define the local energy density 



Ul = ik[ + 2^177 + - xi+i) 2 ■ Using the current con- 
servation equation du/dt + dj/dx = and the equations 
of motion we then find that ji — (xixi-i + xi-\Xi)/2. 
The steady state current can now be computed by using 
the explicit solution in Eq. (7). We get 

/oo 
-oo 

+k' 2 Z^)Z l -\ N (-cj)I'(u;)} (8) 

The matrix Z is tridiagonal and using some of its special 
properties [9] we can reduce the current expression to the 
following simple form: 



ft) - ^ / 
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(/-/') 
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where J(uj) = k 2 k' 2 u)b(u))b' '(uj) / V | Y^jv | 2 has the phys- 
ical interpretation as the total heat current in the wire 
due to all right-moving (or left-moving) scattering states 
of the full Hamiltonian (system + reservoirs). Such scat- 
tering states can be obtained by evolving initial unper- 
turbed states of the reservoirs with the full Hamiltonian 
[9,10]. We have denoted by Yj ;TO the determinant of the 
submatrix of Z beginning with the Ith row and column 
and ending with the mth row and column. Similarly let 
Di, m denote the determinant of the submatrix formed 
from <F 

For the special case when the reservoirs are also one di- 
mensional chains with nearest neighbor spring constants 
Kim = 1 an d the coupling constants k, k' are set to unity, 
we have: G+ p = G plp , = e~' lk where ui = 2sin(fc/2), 
I(uj) = /(w)sin(fc)/7r for | u |< 2 and I(uj) = for 
| co |> 2. In this case Eq. (9) simplifies further and has 
an interpretation in terms of transmission coefficients of 
plane waves across the disordered system. We get: 



J=^y du(j\t N (w)\ 2 (f - /') where (10) 
, 2 4sin 2 (fc) 



e ik {D 2 , Ns + £>i,jv.-i) + e l2k D 2 .N s -i\ 



is the transmission coefficient at frequency ui. We have 
thus obtained the Landauer formula [11] for phononic 
transport. It is only in this special case of a one- 
dimensional reservoir and perfect contacts that we get 
the Landauer formula. The reason is that, only in this 
case is the transmission through the contacts perfect, 
and this requirement is one of the crucial assumptions 
in the Landauer derivation. Note that in Eq. (10) (i)the 
transmission coefficient does not depend on bath prop- 
erties and (ii) transmission is only through propagat- 
ing modes. For general reservoirs where we need to use 
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Eq. (9) the factor J(ui) involves not just the properties 
of the wire but also the details of the spectral functions 
of the reservoirs. Thus the conductivity of a sample can 
show rather remarkable dependence on reservoir prop- 
erties as we shall see below. The above Landauer-like- 
formula has earlier been stated in [12] and derived more 
systematically in [13]. We note that in the high temper- 
ature limit T , T" — > oo Eq. (10) reduces to the classical 
limit obtained exactly in [1,5,6]. 

Asymptotic system-size dependences : In the case of 
electrical conduction the conductance of a long disor- 
dered chain decays exponentially with system size as a 
result of localization of states. In the case of phonons the 
long wavelength modes are not localized and can carry 
current. This leads to power-law dependences of the cur- 
rent on system size as has been found earlier in the con- 
text of heat conduction in classical oscillator chains. A 
surprising result is that the conductivity of such disor- 
dered chains depend not just on the properties of the 
chain itself but also on those of the reservoirs to which 
it is connected. It can be shown [1] that the asymptotic 
properties of the integral in Eq. (9) depend on the low 
frequency (w < 1/N 1 / 2 ) properties of the integrand. This 
means that we will get the same behaviour as in the clas- 
sical case. The details can be found in [1,5,6]. Here we 
summarize some of the main results: 

(i) For the one-dimensional reservoir, which corre- 
sponds to the classical Rubin-Greer model [6], we get 

J ~ ]vT72- 

(ii) The case of reservoirs which give delta-correlated 
Langevin noise correspond to the classical Lebowitz 
model [14,5] where one gets J ~ jjkj2- 

(hi) In general one gets J ~ where a depends 
on the low frequency behaviour of the spectral functions 
Gpp(w) and Gy p /(u/) [1]. 

Note that the case a < 1 leads to infinite thermal 
conductivity while a > 1 gives a vanishing conductivity. 
Thus, depending on the properties of the heat baths, the 
same wire can show either a vanishing thermal conductiv- 
ity or a diverging one. The usual Fourier's law would pre- 
dict J ~ 1/N, independent of reservoirs. Thus Fourier's 
law is not valid in quantum harmonic chains, even in the 
presence of disorder. This breakdown of Fourier's law in 
ID systems has been noted in a number of earlier studies 
on classical systems [4] which have looked at the effects 
of scattering both due to impurities and nonlincarities. 

Temperature profiles: The local temperature of a par- 
ticle can be determined from its average kinetic energy , 
kei = (pf/(2mi)). We get 

1 f x 

kei = - dLomiLo 2 [k 2 Z-^{co)Z^(-Lu)I(Lu) 

2 J-oo 

+k l2 Z-}(u J )Z-}(-u J )I l {uj)] (11) 

This is straightforward to evaluate numerically for given 
systems and reservoirs. Here we shall consider only the 



special case of heat transmission through a perfect one- 
dimensional harmonic chain attached to one-dimensional 
reservoirs. For the case of perfect contacts (i.e k = k' = 
1) Eq. (11) simplifies (for large N) to 



kei 



1 

8^ 



dkuj k [f(uj k ) + /'(wfc)]. 



(12) 



For T = T' one can verify that this corresponds to the 
equilibrium kinetic energy density on an infinite chain. 
On the other hand if we make the couplings weak (i.e. 
k, k' << 1 ) then we can verify that the kinetic energy 
density profile for T = T' corresponds to that of a fi- 
nite chain of length TV in thermal equilibrium ( in the 
quantum case the equilibrium energy density depends on 
I). 

In the non equilibrium case, at high temperatures, the 
local temperature is given by 7} = 2kei and from Eq. (12) 
we get Ti = (T + T')/2 which is the classical result [15]. 
At low temperatures and imperfect contacts k, k' ^ 1 
we evaluate the local kinetic energy profile numerically 
using Eq. (11). As can be seen in Fig. 1 the temperature 
in the bulk still has the same constant value. At the 
boundaries however we see a curious feature noted earlier 
by [15,7]: the temperature close to the hot end is lower 
than the average temperature while that at the colder 
end is higher than the average. The temperature profiles 
at lower temperatures are only quantitatively different. 



Imperfect contacts 
Perfect contact 



FIG. 1. Kinetic energy density profile in a pure harmonic 
wire (N=64) attached to one dimensional reservoirs at tem- 
peratures T = 1 (left) and T' = 0.5 (right), for perfect and 
imperfect contacts (k = k' = 0.9). The temperatures consid- 
ered are quite low and so the bulk temperature is different 
from the classically expected value T av = 0.75. 

In summary we have used the formulation of Ford- 
Kac-Mazur to study heat transport in disordered quan- 
tum harmonic chains. Exact formal relations for the heat 
current have been derived and it is shown that one ob- 
tains Landauer-like results for special reservoirs. Earlier 
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results on classical chains are obtained as limiting cases. 
Finally we make a couple of predictions that are interest- 
ing from the experimental point of view: (i) the large sys- 
tem size behaviour of the heat current is a power law and 
the power depends on reservoir properties (ii) tempera- 
ture profiles in perfect wires show somewhat counterintu- 
itive features close to contacts. A number of experiments 
on heat transport in insulating nanowires have recently 
been carried out [3]. It would be interesting to see if our 
predictions, which are true for strictly one-dimensional 
chains, can be verified in such experiments. 

I thank Ulrich Zurcher for pointing out errors in an 
earlier version of this paper. 
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